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ABSTRACT 


Recognizing  the  uncertainties  as  to  future  available  fuel  sources  and  char¬ 
acteristics,  A.F.O.S.R.  established  a  Center  of  Excellence  in  Combustion  Related 
Fuel  Problems.  The  program  concentrates  on  chemical  aspects  of  combustion  prob¬ 
lems.  This  document  is  the  Annual  Report  for  the  period  1  October  1979  to  30 
September  1980  and  reports  solely  on  the  research  progress  made. 

Results  on  a  detailed  and  overall  reaction  mechanism  for  propane  pyrolysis 
have  been  published.  Some  of  the  first  results  on  the  high  temperature  oxida¬ 
tion  of  aromatics  have  been  obtained.  Most  of  this  work  has  been  with  toluene 
and  a  general  mechanism  of  toluene  oxidation  has  been  proposed.  The  important 
intermediate  formed  is  benzaldehyde  and  the  process  evolves  into  the  oxidation 
of  formyl  and  phenyl  radical.  The  indications  are  that  all  the  alkylated  aro¬ 
matics  follow  the  sequence  of  initial  oxidation  of  the  akyl  group  and  then  the 
oxidation  of  phenyl  radical.  If  so,  complexities  with  respect  to  overall 
kinetic  modelling  for  aromatics  may  not  exist.  The  excellent  data  obtained  are 
due  in  great  part  to  the  development  of  a  high  temperature  sampling  system  for 
the  Princeton  flow  reactor  and  this  system  is  a  significant  contribution  in  it¬ 
self. 

Extensive  results  are  reported  on  the  sooting  tendencies  of  fuels.  In  par¬ 
ticular  the  effects  of  oxygen  index  and  temperature  on  sooting  diffusion  flames 
have  been  evaluated.  Results  verify  our  earlier  postulate  that  it  is  the  stoi¬ 
chiometric  flame  temperature  and  the  fuel  structure  that  determine  sooting  ten¬ 
dencies  in  diffusion  flames.  The  temperature  based  on  the  overall  stoichiometry 
based  the  flow  rates  of  fuel  and  oxidizer  is  not  important.  The  oxygen  concen¬ 
tration  in  the  far  field  as  it  affects  burn-up  rate  is  important. 

Combustion  phenomena  inherently  involves  a  complex  interaction  of  many  sim¬ 
ultaneous  reactions  and  associated  transport  phenomena.  Efforts  to  develop  a 


powerful  family  of  sensitivity  analysis  tools  to  quantitatively  address  the  com¬ 
plex  interaction  between  a  combustion  system's  chemical  and  transport  components 
are  detailed. 


> 


I .  INTRODUCTION 


The  importance  of  chemically-related  combustion  problems  has  led  to  the  es¬ 
tablishment  of  an  integrated  fuels  research  program  under  A.F.O.S.R.  support. 

This  document  is  the  annual  report  of  progress  under  that  effort  for  the  contract 
year  1  October  1979  to  30  September  1980. 

Research  emphasis  is  entered  on  the  pyorlysis  and  fuel -rich  oxidation  of  the 
aliphatic  and  aromatic  hydrocarbons,  the  controlling  mechanisms  of  soot  formation 
and  destruction,  and  sensitivity  analysis  techniques  as  applied  to  combustion 
problems.  The  progress  made  and  results  obtained  during  the  subject  contract 
year  are  reported  in  the  following  sections. 

II.  HYDROCARBON  OXIDATION  AND  PYROLYSIS 

Work  on  the  pyrolysis  of  propane  was  completed  before  the  contract  period 
began,  however  the  writing  and  publication  of  the  material  proceeded  into  the 
period  of  concern.  The  paper  entitled  "An  Overall  and  Detailed  Kinetic  Study  of 
the  Pyrolysis  of  Propane"  by  D.J.  Hautman ,  R.J.  Santoro,  F.L.  Dryer  and  I.  Glassman 
was  accepted  for  publication  in  the  International  Journal  of  Chemical  Kinetics1. 

After  the  completion  of  this  work,  our  efforts  concentrated  on  the  very  im¬ 
portant  problem  of  the  oxidation  of  aromatics.  OSHA  regulations  at  the  time  the 
study  began  forbade  the  use  of  benzene,  consequently  our  initial  efforts  concen¬ 
trated  on  toluene. 

As  all  previous  oxidation  studies,  the  work  was  performed  on  the  Princeton 
turbulent  flow  reactor.  The  use  of  the  aromatics  in  the  reactor  posed  substan¬ 
tial  problems  in  sampling  and  analysis  of  species  formed  during  the  oxidation  pro¬ 
cess.  Whereas,  in  the  oxidation  of  aliphatics,  experience  had  shown  that  the 
species  formed  were  low  boiling  point  compounds  whose  peaks  in  gas  chromatographic 
analysis  were  readily  known;  it  was  to  be  expected  in  the  aromatic  work  that  much 


< 


2 


heavier  species  would  be  found. 

The  presence  of  unknown  high  boiling  point  species  necessitated  the  develop¬ 
ment  of  a  unique  sample  storage  technique.  The  development  of  this  technique 
and  the  successful  functioning  and  use  of  our  GC/MS  apparatus  were  major  accom¬ 
plishments  during  the  subject  year,  but  more  importantly,  permitted  us  to  obtain 
sufficient  data  on  the  oxidation  of  toluene  to  begin  to  understand  its  complex 
mechanism. 

A.  New  Sampling  System 

Both  temperature  and  species  profiles  are  measured  along  the  length  of  the 

reactor  duct.  The  temperature  is  measured,  as  in  earlier  flow  reactor  studies, 

with  a  silica  coated  Pt  30%  Rh/Pt  6%  Rh  alloy  thermocouple  with  a  bead  diameter 

2 

of  .0015"  and  a  spatial  resolution  of  about  twenty  bead  diameters.  Dryer  has 
discussed  this  technique  at  length  in  an  early  Ph.D.  thesis  under  A.F.O.S.R. 
sponsorship . 

The  gas  sampling  techniques,  in  principle,  are  also  those  used  in  earlier 
flow  reactor  studies:  combustion  gases  are  withdrawn  from  the  reactor  with  a 
water  cooled  probe  and  are  stored  for  later  analysis  by  a  gas  chromatograph.  A 
major  purpose  of  the  present  work,  however,  was  to  develop  a  sample  acquisition 
and  storage  system  which  operates  at  temperatures  above  ambient  to  prevent  con¬ 
densation  and  polymerization  reactions.  This  system,  which  replaces  the  ambient 
temperature  glass  bottle  storage  system  previously  in  use,  consists  primarily  of 

1.  a  variable  temperature  probe  coolant  system  for  quenching  of  the  combustion 
gases  at  a  controlled  temperature 

2.  a  heated  multiposition  sample  storage  valve  with  sixteen  12  cc .  storage  loops 
heated  gas  transfer  equipment  which  comprises  the  interface  between  the  mul¬ 
tiposition  valve  and  either  the  flow  reactor  or  a  gas  chromatograph 
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The  final  design  of  the  heated  sampling  system  was  the  result  of  essentially 
two  stages  of  development.  The  first  stage  focussed  on  the  problems  associated 
with  high  temperature  acquisition  and  storage  of  samples  and  resulted  in  the 
choice  of  a  multiposition  valve  as  the  principal  storage  device.  As  development 
of  this  first  generation  system  proceeded,  however,  several  problems  were  iden¬ 
tified,  all  of  them  associated  with  operational  characteristics  of  the  valve.  The 
second  generation  system,  through  a  more  indirect  approach  to  the  storage  problem, 
circumvented  the  valve  limitations  and  retained  the  essential  features  of  the 
first  design.  The  fundamental  problems  of  high  temperature  storage  and  the  op¬ 
eration-related  valve  problems  were  all  minimized  or  eliminated  in  the  second  sys¬ 
tem.  Both  systems  are  discussed  below  because  an  understanding  of  the  system's 
evolution  will  provide  insights  into  the  details  of  its  final  design  and  into  the 
limitations  of  its  use. 

Successful  chemical  sampling  involves  satisfactory  accomplishment  of  each  of 
the  following  goals: 

1.  isolation  of  a  flow  of  reactive  media  from  a  specific  point  in  the  reaction 
environment 

2.  prevention  of  further  chemical  or  thermodynamic  alteration  of  the  contained 
species 


introduction  of  a  defined  quantity  of  this  sample,  intact,  to  an  analytical 


instrument 


The  room  temperature  system  previously  used  in  our  A.F.O.S.R.  kinetic  studies 


reasonably  satisfied  these  requirements  (for  known  stable  species  except  form¬ 
aldehyde).  It  incorporated  storage  of  samples  in  glass  bottles  at  room  tempera¬ 
ture  and  required  pressurization  of  the  samples  (for  injection  onto  a  gas  chrom¬ 
atograph)  by  displacement  of  a  sample  volume  with  mercury. 
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Condensation,  which  was  not  a  problem  in  earlier  flow  reactor  studies 
(because  the  fuels  considered  in  those  studies  gave  rise  to  intermediate  species 
and  products  which  were  gases  or  high  vapor  pressure  liquids  at  room  temperature)  , 
was  confirmed  early  as  a  problem  in  the  aromatics  studies.  First,  a  yellow-brown 
gummy  film  was  noted  on  the  outer  surface  of  the  sampling  probe  following  aro¬ 
matic  oxidation  runs.  Second,  a  carbon  balance  through  the  reaction  for  aromatic 
oxidation  data  gathered  with  the  older  sampling  system  showed  a  loss  in  total 
carbon  content  in  the  intermediate  stages  of  the  reaction,  where  high  boiling 
point  intermediates  are  likeliest  to  exist.  And  third,  an  examination  of  the 
available  literature  revealed  postulated  and  observed  intermediates  of  critical 
importance  which  have  extremely  small  vapor  pressures  at  the  quench  temperatures 
used  with  the  older  sampling  system.  Though  it  is  impossible,  a  priori,  to  theo¬ 
retically  estimate  the  magnitude  of  a  condensation  problem  --  the  intermediates 
and  their  concentrations  would  need  to  be  known  for  such  a  calculation  --  the 
heated  sampling  and  storage  approach  suggested  here  seemed  reasonable  for  two 
reasons.  First,  it  allows  for  positive  determination  of  whether  or  not  a  conden¬ 
sation  problem  exists;  and  second,  it  is  designed  for  easy  upgrading  to  higher 
temperatures  if  a  problem  is  identified  and  persists  with  the  present  system. 

The  system  described  here,  at  the  very  least,  allows  for  definitive  analysis 
of  the  condensation  problem.  This  was  not  possible  with  the  room  temperature  sys¬ 
tem  because  the  condensation  limited  partial  pressures  at  ambient  temperatures  of 
some  possible  critical  intermediates  (e.g.  hydroquinone  and  pyrocatechol)  are 
below  the  detection  limits  of  the  flame  ionization  detector  (FID)  used  for  the 
analysis.  In  other  words,  at  room  temperature,  insufficient  amounts  of  these  in¬ 
termediates  remain  in  the  gas  phase  for  detection  (for  sample  sizes  appropriate 
to  gas  chromatographic  analysis  of  other  species).  At  the  higher  quench  tempera¬ 
tures  of  the  heated  system,  this  is  not  the  case,  and  these  high  boiling  point 
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species  can  be  identified  and  quantified.  Comparison  of  the  amount  detected  with 
the  vapor  pressure  of  the  species  at  the  sampling  temperature  will,  in  turn,  allow 
determination  of  the  species,  if  any,  which  are  condensing.  The  system  can  then 
be  upgraded  to  a  point  where  the  tradeoff  between  the  condensation  of  stored 
gases  and  their  reactivity  is  optimized. 

There  are  two  advantages  to  controlling  condensation  in  this  way.  First, 
major  high  boiling  point  intermediate  species  can  be  accurately  measured.  But 
additionally,  minor  aromatic  species  —  which  will  not  greatly  affect  the  carbon 
balance  or  an  "overall"  analysis  of  the  data  --  can  be  detected,  identified,  and 
quantified.  The  heated  sampling  system  interfaces  with  our  GC/MS  in  a  manner 
completely  analagous  to  its  GC  interface,  and  allows  ready  identification  of 
species  which  are  impractical  to  identify  with  a  GC  alone  and  which  would  be  lost 
to  condensation  at  room  temperature. 

The  storage  temperature  of  70°C  has  a  second  advantage  related  to  the  in¬ 
tegrity  of  the  stored  sample:  by  preventing  polymerization  reactions  associated 
with  the  storage  of  formaldehyde  at  room  temperature,  it  should  allow  accurate 
formaldehyde  analysis. 

Also,  the  heated  sampling  system  has  the  advantages  of  preventing  exposure 
of  stored  samples  to  light  and  to  adsorptive  surfaces  (e.g.  mercury,  systhetic 
rubber)  ,  both  of  which  occurred  in  the  older  system.  All  surfaces  which  contact 
sample  gases  in  the  heated  sampling  system  are  stainless  steel,  so  the  exposure 
to  light  that  resulted  from  storage  in  glass  bottles  on  an  exposed  sampling  rack 
is  circumvented.  Furthermore,  gases  are  pressurized  for  injection  with  a  stain¬ 
less  steel  welded  bellows  pump  and  not  with  the  mercury  system  previously  in  use. 

The  heating  of  a  number  of  stored  samples  warrants  a  more  compact  system 
than  the  glass  bottles  previously  in  use.  Further,  the  transfer  and  handling  of 
individual  hot  samples  for  analytical  purposed  requires  either  the  use  of  a  sep- 


arate  heated  and  insulated  oven  for  each  sample  or  a  method  of  remote  control  of 
a  collection  of  samples  in  a  single  heated  environment.  A  remotely  actuated  mul¬ 
tiposition  valve  with  several  storage  loops  satisfies  both  of  these  requirements 
(Figure  1) . 

The  most  direct  approach  to  the  use  of  such  a  valve  as  the  storage  unit  for 
a  gas  sampling  requires  a  system  which  consists  primarily  of  three  components: 
the  probe  coolant  system,  the  gas  transfer  equipment  (including  the  lines,  the 
pump  and  the  manostat) ,  and  the  sample  storage  system  (the  multiposition  valve) . 
Wherever  necessary,  the  components  are  heated  to  prevent  condensation. 

Our  sampling  probes  are  constructed  of  304  stainless  steel  and  are  cooled  by 
two  water  jackets,  one  inside  the  other.  The  temperature  of  the  water  passing 
through  the  jackets  at  the  probe  inlet  can  be  varied  between  about  15°C  and  80°C. 
Through  control  of  this  water  temperature,  the  temperature  to  which  the  combustion 
gases  are  quenched  can  be  controlled. 

The  probe  coolant  system  consists  simply  of  a  recirculating  water  stream 
which  is  successively  cooled  by  a  heat  exchanger  and  heated  by  the  hot  reactor 
(through  the  probe) .  Since  heat  addition  to  the  coolant  stream  varies  as  the 
probe  is  extended  into  the  reactor,  variable  heat  transfer  at  the  heat  exchanger 
is  required  to  produce  the  isothermal  quench  temperature  chosen  for  a  study.  This 
is  accomplished  with  a  remotely  actuated  valve  that  varies  the  secondary  flow  rate 
through  the  heat  exchanger. 

Operation  of  the  first  generation  sampling  system  was  as  follows:  combustion 
gases  were  withdrawn  from  the  reactor  through  the  probe  and  a  length  of  heated 
line  by  a  heated  stainless  steel  welded  bellows  pump.  The  bellows  pump  raised 
the  pressure  of  the  sampled  gases  above  atmospheric  pressure,  and  the  flow  pro¬ 
ceeded  through  a  heated  7  micron  filter  and  a  temperature  equilibration  coil  to 
the  sample  storage  valve.  The  temperature  equilibration  coil  and  the  sixteen  loop 


storage  valve  were  contained  in  a  separate,  portable,  thermostatted  oven  assembly. 
The  pressure  of  the  gases  being  sampled  was  controlled  by  a  manostat /vacuum  oump 
system  downstream  of  the  valve.  Each  of  the  valve  loops  had  a  pre-determined 
volume  of  .5  cc .  When  the  valve  was  actuated  and  the  flow  switched  from  one 
loop  to  the  next,  a  defined  quantity  of  sample  was  trapped  in  the  first  loop 
(P,  V  and  T  are  all  defined).  Up  to  fifteen  samples  could  be  stored  in  this  way, 
with  the  remaining  loop  open  to  the  flow  path.  Once  filled,  the  multiposition 
valve  system  was  carried  (in  its  oven)  to  a  gas  chromatograph,  where  the  loops 
were  flushed,  one  by  one,  into  the  GC  columns  for  analysis.  In  this  design,  the 
pressure,  volume  and  temperature  of  the  stored  sample  were  chosen  to  accommodate 
gas  chromatographic  sample  size  requirements. 

Thissystem  fulfilled  in  principle  each  of  the  sampling  system  requirements 
outlined:  defined  quantities  of  sample  were  isolated  and  trapped  in  the  multi¬ 
position  valve  loops;  the  storage  temperature  was  low  enough  to  prevent  signifi¬ 
cant  reaction,  but  high  enough  to  prevent  condensation  or  polymerization ;  and 
the  samples  were  transferred  directly  to  a  gas  chromatograph  for  analysis.  No¬ 
where  in  the  transfer  or  storage  of  samples  were  they  exposed  to  light  or  made  to 
contact  adsorptive  substances  like  mercury  or  synthetic  rubber.  (The  adsorptive 
sites  on  the  stainless  steel  in  the  system  were  most  likely  pacified  during  the 
flushing  of  sample  through  the  system  during  sample  acquisition).  Unfortunately, 
operational  limits  of  the  multiposition  valve  give  rise  to  several  problems  out¬ 
lined  below: 

1.  Loop  to  Loop  Volume  Variations  -  The  volumes  of  the  sixteen  loops  were  dif¬ 
ficult  to  match,  so  sample  sizes  vary  slightly.  Calibrating  each  loop  in¬ 
dividually  could  have  eliminated  this  problem,  but  at  the  cost  of  additional 


8 


il 

2.  Leakage  into  or  out  of  the  Loops  or  Valve  Body  -  No  leaks  through  individual 

loops  were  identified,  but  a  small,  progressive  loss  of  sample  over  time  was 
observed  when  samples  were  stored  at  above  atmospheric  pressure.  Leakage  I 

out  made  determination  of  exact  sample  size  impossible.  The  inverse  prob-  i 

lem  of  leakage  in  for  storage  at  sub-atmospheric  pressure  lead  to  dilution 

and  distortion  of  particular  species  (especially  oxygen). 

3.  Destruction  of  Stored  Samples  due  to  Passage  of  High  Pressure  GC  Carrier  Gas 
Through  the  Multiposition  Valve  -  This  was  a  large  effect.  Since  the  GC  car¬ 
rier  gas  must  enter  the  chromatographic  columns  at  about  60  psig  in  order 

j 

■1 

to  ;  rcome  the  resistance  of  the  columns,  helium  at  60  psig  must  pass  through 
the  multiposition  valve  in  order  to  flush  the  samples  onto  the  columns.  Al-  ? 

though  the  manufacturer  specifies  that  the  valve  sealing  is  good  to  400  psig,  j 

we  found  that  flushing  the  samples  with  high  pressure  helium  caused  signifi-  : 

cant  distortion  of  four  out  of  fifteen  stored  samples.  The  four  that  were 
destroyed  were  always  the  first  and  last  three  loops  analyzed,  no  matter  which 
loop  was  chosen  as  the  starting  point.  Different  loops  were  chosen  as  the 
starting  point.  The  pattern  and  magnitude  of  the  distortion  was  clear,  though 
the  mechanism  was  uncertain. 

4.  Diffusion  of  Sample  Between  Adjacent  Loops  -  The  manufacturer  of  the  multi- 
position  valve  specified  a  diffusion  rate  of  1.0  x  lO’^torr-liters/sec .  of 
helium,  a  value  roughly  substantiated  by  tests  conducted  in  this  laboratory. 

For  the  sample  sizes  chosen  for  accurate  chromatographic  analysis  (about  .5 
cc.  -  atm.),  and  for  storage  times  (25+  hours)  required  for  analysis  of  fif¬ 
teen  samples  containing  the  species  expected  in  aromatics  oxidation  reactions,  t 

this  small  rate  could  significantly  alter  stored  samples.  The  diffusion  dis-  j 

tortion  was  demonstrated  by  filling  the  sample  valve  loops  with  either  cali¬ 
bration  gas  (CH^,  C ,  C2H(3»  C^Hg,  C'0,  or  Pure  N->  anc*  analyzing 
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the  stored  samples  over  time.  Errors  in  excess  of  15%  were  measured  in  less 
than  25  hours. 

These  four  problems,  especially  the  last  three,  made  this  critical  approach 
to  storage  and  analysis  of  samples  unacceptable.  The  design  philosophy  had  to  be 
changed,  and  in  the  second  generation  system,  the  errors  resulting  from  the  prob¬ 
lems  discussed  above  were  either  eliminated  altogether  or  made  negligibly  small. 

Figure  2  shows  a  schematic  of  the  new  system.  Again,  samples  are  drawn  through 
the  probe  with  a  heated  bellows  pump  and  stored  in  loops  on  a  sixteen-position  valve. 
But  the  volume  of  each  of  the  loops  has  been  increased  from  .5  cc.  to  12  cc.  (Fig¬ 
ure  3),  and  the  storage  pressure  is  now  about  2  atmospheres  instead  of  .5  atmos¬ 
pheres.  The  higher  pressure  is  achieved  by  simply  replacing  the  manostat/vacuum 
system  of  the  previous  design  with  a  needle  valve,  which  throttles  the  flow  to 
achieve  the  desired  value.  Although  throttling  causes  a  drop  in  the  flow  rate 
through  the  system,  less  than  two  sweep  volumes  are  required  to  ensure  complete 
displacement  of  any  residual  gases  in  the  transfer  lines,  and  the  system  is  de¬ 
signed  to  achieve  this  easily.  Both  flow  rate  and  pressure  are  measured  at  the 
outlet  of  the  multiposition  valve. 

The  new  system  stores  almost  50  times  the  sample  size  appropriate  for  GC 
analysis.  A  less  direct  approach  to  injection  of  samples  onto  the  chromatographic 
columns  was  therefore  necessary.  A  secondary  1  cc.  sample  loop  is  now  used  to 
set  the  sample  size  for  analysis.  This  sample  loop  is  connected  on  the  inlet 
side  to  the  multiposition  valve  (i.e.  to  the  stored  samples)  through  a  length  of 
heated  line,  and  on  the  outlet  side  to  a  manostat/vacuum  system.  When  the  multi¬ 
position  valve  is  actuated  and  a  new  storage  loop  is  exposed  to  the  sample  line, 
the  gases  in  the  12  cc.  loop,  which  are  stored  at  2  atmospheres,  flow  through  the 
1  cc .  secondary  sample  loop  until  the  pressure  in  the  whole  system  is  equilibrated 
to  .5  atmospheres  by  the  manostat.  Sufficient  excess  sample  is  stored  in  the  mul- 
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tiposition  valve  loops  to  assure  complete  displacement  of  any  residual  gases  in 
the  secondary  sample  valve  and  transfer  lines.  The  .5  cc.-atm.  sample  in  the 
secondary  sample  valve  is  then  injected  onto  the  gas  chromatograph  columns  for 
analysis.  Actuation  of  both  the  multiposition  valve  and  the  secondary  sample 
valve  is  controlled  by  a  microprocessor  indigenous  to  the  gas  chromatograph,  so 
fully  automatic  analysis  of  successive  samples  is  possible.  This  procedure  not 
only  assures  a  constant  systematic  error,  but  minimizes  sample  storage  times 
and  reduces  dramatically  the  manpower  which  must  be  devoted  to  the  analysis  of 
samples.  The  previous  (manual)  injection  technique  required  operator  attention 
throughout  the  25+  hour  analysis  time.  Once  the  present  system  is  interfaced 
with  the  gas  chromatograph  and  the  GC  microprocessor  is  programmed,  nj  additional 
operator  attention  is  required  for  complete  analysis. 

The  importance  modification  introduced  in  this  second  generation  design  is 
the  separation  of  the  sample  storage  requirements  from  the  GC  sample  size  require¬ 
ments.  This  separation  allows  a  flexibility  in  parameters  that  virtually  elimi¬ 
nates  the  problems  previously  discussed. 

Because  the  sample  injected  onto  the  GC  columns  is  set  in  the  same  loop  each 
time,  the  effects  of  loop  volume  variability  are  eliminated.  The  sample  leakage 
problems  are  also  circumvented:  since  samples  are  stored  at  a  positive  pressure, 
gases  leak  out,  not  in.  Loss  of  sample  at  the  small  rates  indicated  by  earlier 
tests  results  in  loss  of  a  fraction  of  the  ste^d  sample.  This  loss  is  not  re¬ 
flected  in  the  sample  set  in  the  secondary  sampling  valve,  since  a  large  excess 
of  sample  gas  remains  in  the  storage  loops  under  all  conditions.  Any  difficulties 
associated  with  passage  of  high  pressure  GC  carrier  gas  through  the  multiposition 
valve  are  also  eliminated.  Finally,  the  diffusion  problem  is  reduced  to  a  negli¬ 
gible  level.  Although  the  absolute  diffusion  rate  remains  the  same,  the  much 
larger  quantity  of  sample  stored  in  each  loop  dilutes  the  contamination  by  a  fac- 
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tor  of  about  fifty. 

Several  tests  were  conducted  to  ascertain  the  magnitude  of  the  persisting 
diffusion  error.  In  these  tests,  as  in  earlier  checks  on  diffusion  errors,  loop.' 
were  alternately  filled,  at  pressures  expected  during  routine  operation,  with  cal¬ 
ibration  gas  and  with  pure  nitrogen.  The  samples  were  then  analysed  in  exactly 
the  way  they  would  be  during  actual  operation.  The  diffusion  of  calibration  gas 
into  the  pure  nitrogen  loops  was  then  plotted  as  a  function  of  time.  The  dif¬ 
fusion  rate  was  calculated  to  be  about  3x10  ^  torr-liters  He/sec.,  in  rough  agree¬ 
ment  with  the  manufacturer's  specified  value,  but  the  distortion  of  the  stored 
sample  is  shown  to  be  less  than  1%,  even  for  the  lightest  gas  expected  in  our 
hydrocarbon  oxidation  studies  (CH^) .  Note  that  since  diffusion  rates  are  inversely 
proportional  to  the  square  root  of  the  molecular  weight  of  the  diffusing  species, 
errors  for  most  species  will  be  significantly  smaller. 

In  summary,  the  advantages  of  the  sampling  system  presented  here  are: 

1.  control  of  the  sample  acquisition  and  storage  temperatures  (up  to  about  S0C 

C) 

2.  sample  isolation  from  adsorptive  surfaces 

3.  sample  isolation  from  light  sources 

4.  automation 

B.  Toluene  Oxidation  Results 

Species  profiles  typical  of  those  obtained  in  these  early  toluene  oxidation 
experiments  are  given  in  Figures  4  and  5.  All  profiles  are  given  in  the  thesis 
by  Euchner^.  Two  plots  are  given  for  each  reaction  studied.  Note  that  the  mole 
fraction  of  carbon  is  plotted  in  these  figures,  not  the  species  mole  fraction. 

Seven  times  the  mole  fraction  of  toluene  (C^Hg)  is  plotted,  for  example,  while  the 
mole  fraction  of  CO  is  plotted  directly.  The  equivalence  ratio  for  these  experi - 


ments  was  varied  between  about  .Ob  and  .43,  and  the  temperature  ranged  between 
10SD  k  and  1150  k.  The  major  intermediate  species  were  identified  and  quantified, 
t  iough  the  quantification  of  minor  aromatic  species  remains  uncertain.  These  f.t' 

MS  results  confirmed  the  identifications  performed  with  the  GC  and  gave  additional 
information  on  species  that  could  not  be  chromatographical ly  detected.  The  GC 
data,  which  report  relative  mole  fractions  of  the  major  intermediates  and  peak 
reaction  coordinate  data  for  both  the  major  and  minor  species,  give  some  indica¬ 
tion  of  the  mechanism.  Note  first  the  general  sequential  and  overlapping  trend: 

TOLUENE  -  AROMATIC  -  ALIPHATIC  -  CO  CO., 

INTERMEDIATES  INTERMEDIATES 


Within  the  aromatic  intermediates,  benza Idchyde  and  benzene  are  the  dominant 
species,  and  phenol,  styrene  ( or  the  xylenes  l,  ethylbenzene  and  benzyl  alcohol 
are  formed  in  lesser  amounts. 

The  data  indicate  that,  for  all  conditions  studied,  the  major  courses  of  tin 
reaction  is  through  benza  ldehyde  to  benzene,  followed  by  ring  fragmentation.  Hu 
most  probable  mechanism  based  on  this  observation  is  given  below.  (R  is  am  radi¬ 
cal  --  a  hydrocarhon  radical,  H,  C ,  CH  or  am  other  radical'. 


(1) 


(2) 


(3) 


INITIATION 


CHAIN 

PROPOGATION 


CHAIN 
BRANCH  INC 


(4) 


CHO 


Iinti.it  ion  occurs  primarily  by  abstraction  of  an  alkyl  hydrogen  from  toluene 
1  he  strength  of  a  l.'-H  bond  in  the  methyl  group  is  about  S5  kcal/mole,  significan- 
t 1>  less  than  the  strength  of  the  bond  between  the  methyl  group  and  the  ring  (102 
kcal/mole!  or  the  strength  of  a  ring  C-H  bond  (112  Real /mole j. 

Most  fuel  molecules  probably  then  go  through  some  variation  of  the  chain  pro 
pagating  step.  Reaction  (2),  which  leads  to  a  large  benzyl  radical  pool. 

Chain  branching  occurs  through  Reaction  (3).  Though  endothermic,  this  reac¬ 
tion  appears  to  be  the  only  feasible  one  at  this  temperature. 

Unimolecular  decomposition  of  the  benzoxy  radical  to  benzaldehyde  is  postu¬ 
lated  in  Reaction  (4),  a  step  suggested  by  the  large  quantities  of  benzaldehyde 
and  small  quantities  of  benzyl  alcohol  detected  in  this  study. 

The  benzaldehyde  formed  in  Reaction  (4)  could  react  as  in  Reactions  (51. 
Radical  attack  is  expected  to  dominate,  especially  at  lower  temperatures.  The 
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aromatic  radical  thus  formed  {  CO  )could  decompose  easily  to  form  the  phenyl 

radical  and  CO.  Alternately,  especially  at  higher  temperatures,  benzaldehyde 
could  interact  with  a  second  body  and  decompose  into  a  phenyl  radical  and  a  formyl 
radical.  The  latter  quickly  decomposes  further  into  a  hydrogen  radical  and  CO, 
as  shown.  Similar  reactions  are  well  established  in  the  oxidation  of  aliphatic 
hydrocarbons.  As  the  temperature  was  increased  in  various  runs,  the  amount  of 
benzene  relative  to  benzaldehyde  is  increased.  This  trend  probably  results  from 
an  increase  in  the  thermal  decomposition  of  benzaldehyde  as  the  temperature  is 
increased* . 

Note  that  the  major  products  of  Reactions  (5)  are  the  phenyl  radical  and  CO. 
The  very  reactive  phenyl  radical  might  easily  abstract  a  hydrogen  atom  to  form 
benzene,  as  shown  in  Reaction  (6).  Note  at  this  point  only  that,  for  six  of  the 
seven  carbon  atoms  in  a  toluene  fuel  molecule  (the  ring  carbon  atoms)  ,  the  oxida¬ 
tion  of  toluene  appears  to  reduce  to  the  oxidation  of  the  phenyl  radical  which 
should  also  occur  in  benzene  oxidation  (though  radical  pools  are  expected  to  be 
different  in  the  two  cases).  Thus  the  important  conclusion  of  the  subject  year's 
work  is  that  the  oxidation  of  toluene  the  methyl  group  is  oxidation  in  a  manner 
similar  to  paraffins  and  the  remaining  phenyl  radical  is  then  oxidized.  Subse¬ 
quently  the  possibility  exists  that  the  understanding  of  the  oxidation  of  the 
alkylated  aromatics  may  hinge  only  on  the  mechanism  of  phenyl  oxidation. 

III.  SOOT  FORMATION  AND  DESTRUCTION 

During  the  previous  contract  year,  we  began  our  initial  experimental  work  on 
sooting  phenomena  by  dealing  with  simple  diffusion  flames  and  measuring  the  so- 


*  Because  of  the  instability  of  benzaldehyde  at  these  temperatures,  further  oxi¬ 
dation  of  benzaldehyde  to  benzoic  acid  is  expected  to  be  unimportant.  It  must  be 
noted,  however,  that  benzoic  acid  would  not  have  been  detectable  in  this  study, 
were  it  to  occur,  because  it  would  have  been  condensed  onto  the  GC  columns  in  the 
analytical  train.  The  fact  that  early  CO.,  (which  would  form  from  the  decomposition 
of  benzoic  acid)  is  not  evident  in  the  flow  reactor  profile  data,  however,  lends 
credence  to  the  supposition  that  little  or  no  benzoic  acid  is  forming. 
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called  smoke-points  or  sooting  heights.  Simultaneous  with  this  effort.  Prof,  ('.la 
man  reviewed  the  field  and  published  an  A.F.O.S.R.  report4  entitled  "Phenomenologi 
cal  Models  of  Soot  Formation  in  Combustion  Systems".  Both  these  efforts  provided 
the  direction  for  the  research  during  the  subject  contract  year.  The  results  ob¬ 
tained  have  been  extensive  and  of  great  practical  consequence. 

Since  initial  efforts  have  now  been  published  only  the  guiding  concepts  will 
be  detailed  here.  Essentially  the  review  led  us  to  the  conclusions  that  in  pre¬ 
mixed  flames  the  dominant  factor  controlling  the  tendency  to  soot  was  the  flame 
temperature,  and  that  the  C/0  equivalence  ratio  in  the  fuel  rich  region  was  only 
a  secondary  factor.  This  conclusion  was  drawn  from  the  early  work  of  Millikin5 
who  postulated  that  the  sooting  tendency  in  pre-mixed  flames  was  a  competition 
between  the  rate  of  formation  of  the  soot  precursors  and  the  rate  of  oxidation 
attack  on  these  precursors.  One  could  show  that  as  the  temperature  was  raised, 
the  oxidative  attack  dominated.  Ke  further  postulated  that  the  trend  for  aromat¬ 
ics  was  different  due  to  their  oxidative  mechanism.  Recent  results^  still  indi¬ 
cate  this  possible,  but  different  from  the  initially  proposed  mechanism. 

For  diffusion  flames  we  postulated  that  the  important  factors  controlling 
the  tendency  of  a  fuel  to  soot  were  the  diffusion  flame  temperature  and  the  fuel 
structure.  However,  in  the  diffusion  flames,  since  there  is  no  oxidative  attack 
on  the  fuel  precursors,  the  higher  the  flame  temperature,  the  higher  the  rate  of 
fuel  pyrolysis  and  precursor  formation,  and  thus  the  greater  tendency  to  soot. 

One  of  our  initial  experimental  efforts  on  sooting  diffusion  flames  was  to 
determine  the  effect  of  inert  additives  mixed  with  the  gaseous  fuels  used.  The 
addition  of  any  inert  additive  lessened  the  tendency  of  the  fuel  to  soot,  but, 
more  importantly,  we  found  that  the  effectiveness  of  any  additive  was  direct ly 
proportional  to  its  molar  specific  heat.  These  results  tended  to  support  the  hy¬ 
pothesis  that  the  temperature  was  the  dominant  factor.  Consequently  we  carried 
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out  very  extensive  controlled  series  of  experiments  during  the  subject  contract 
year . 

In  order  to  place  the  results  obtained  in  the  right  context,  it  is  worth¬ 
while  to  review  the  experimental  procedure  used.  The  data  were  taken  on  a  1  cm 
concentric  fuel  jet  surrounded  by  an  air  stream  flowing  in  a  glass  cylinder  of 
about  10  cm  in  diameter.  The  volumetric  flow  rate  of  the  air  was  sufficiently 
greater  than  the  fuel  flow  rate  so  that  the  system  was  always  in  a  highly  over¬ 
ventilated  mode.  Diffusion  flame  heights  vary  directly  with  the  volumetric  flow 
rate  of  the  fuel.  In  the  concentric  configuration  used  in  these  experiments, 
the  flame  height  varies  slightly  with  the  surrounding  air  flow  until  a  certain 
volumetric  rate  is  reached  and  then  the  flame  height  is  independent  of  the  air 
flow  rate.  Earlier  results  by  other  investigators  are  inconsistent  because  they 
did  not  recognize  f'h:-_,  fc  ;or. 

Every  hydrocarbon  diffusion  flame  appears  luminous  due  to  the  formation  of 
soot.  However  unuer  low  fuel  flow  rates  the  amount  of  soot  formed  is  sufficiently 
small  so  that,  although  the  particles  become  luminous  in  the  high  temperature 
diffusion  flame  reaction  zone,  they  are  burned  as  soon  as  the  flow  streamline 
carry  them  into  the  high  temperature  oxidizer  side  of  the  flame,  as  the  fuel  flow¬ 
rate  increases  the  flame  height  increases.  At  a  given  fuel  flow  rate  sufficient 
soot  apparently  forms  within  the  fuel  jet  so  that  when  the  particles  flow  through 
the  flame  they  cannot  be  oxidized  and  a  streak  of  soot  particles  appears.  In  the 
concentric  configuration  used  an  annulus  opens  up  and  when  the  flame  is  viewed 
frontally,  it  looks  as  if  a  pair  of  wings  appear.  The  flame  height  at  which  this 
phenomenon  occurs  is  called  the  normal  smoke  point  or  sooting  height  and  the 
volumetric  or  mass  flow  rate  corresponding  to  this  height  is  used  as  a  qualitative 
measure  of  the  sooting  tendency.  Thus  the  larger  the  sooting  height  or  volumetric 
or  mass  flow  rate  at  this  height,  the  less  the  tendency  to  soot;  and  vice-versa. 
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The  reason  a  smoke  point  exists  is  at  first  not  clear.  Initially  one  would 
expect  the  larger  the  flame  the  longer  the  pyrolysis  time  and  thus  the  greater 
the  tendency  to  soot.  However  when  the  fuel  jet  size  is  constant,  the  larger 
volumetric  flow  rate  (and  thus  flame  height)  is  obtained  by  increasing  the 
fuel  velocity.  Thus  the  stay  time  of  a  fuel  element  within  the  flame  is  constant, 
no  matter  what  the  height.  Then  the  question  arises  as  to  why  sufficient  soot 
forms  at  certain  height  to  break  through  the  flame.  We  have  concluded  that  as 
the  fuel  flow  rate  is  increased,  more  and  more  air  (though  very  small  quantities) 
is  entrained  from  the  surrounding  stream  near  the  fuel  port  tip  into  the  fuel 
side  of  the  flame.  It  is  well  known  that  small  amounts  of  oxygen  will  catalyze 
the  pyrolysis  of  hydrocarbon  fuels.  Indeed  we  had  observed  this  fact  in  our 
earlier  experiments  with  additives  when  we  doped  our  fuel  with  small  amounts  of 
oxygen.  We  could  reduce  our  smoke  points  appreciably  with  the  addition  of  small 
amounts  of  oxygen;  indeed  the  smoke  point  varied  with  the  amount  of  oxygen  added. 

There  are  important  consequences  of  this  concept.  If  indeed  it  is  the  in- 

4 

duced  oxygen  which  eventually  causes  a  smoke  point,  then  our  earlier  hypothesis 
that  the  pyrolysis  rate  of  the  fuel  at  a  fixed  temperature  is  the  controlling  fac¬ 
tor  in  sooting  tendency  is  correct. 

Experiments  were  then  performed  on  numerous  fuels.  For  each  fuel,  we  deter¬ 
mined  the  smoke  point  as  a  function  of  the  amount  of  inert  diluent  added.  Two 
diluents  were  used  -  nitrogen  and  argon.  Most  of  the  results  were  obtained  with 
nitrogen.  The  results  were  initially  plotted  as  the  volumetric  or  mass  flow  rate 
versus  the  calculated  adiabatic  stoichiometric  flame  temperature  for  the  fuel  with 
the  corresponding  amount  of  diluent  added.  Although  for  diffusion  flames  the 
volumetric  flow  rate  is  the  more  fundamental  one,  the  mass  flow  rate  normalizes 
the  data  for  comparing  fuel  types.  The  plot  of  flow  rate  versus  temperature  gave 
a  series  of  curves  for  each  fuel  that  showed  most  interesting  trends.  There  seemed, 


18 


however,  to  be  a  more  meaningful  manner  to  plot  the  data. 

If  the  fuel  pyrolysis  rate  were  the  controlling  element  in  the  sooting  ten¬ 
dency,  then  the  flow  rate  at  the  smoke  point  (or  the  sooting  height)  was  an  in¬ 
verse  measure  of  the  pyrolysis.  And,  if  the  pyrolysis  rate  followed  Arrhenius 
kinetics,  then  a  plot  of  the  (log  of  the  inverse  of  the  fuel  rate)  versus  1/T 

should  give  straight  lines.  All  the  data  correlated  in  this  manner!  Since  this 

7  8 

work  has  been  published,  ’  we  only  show  the  most  important  graph  in  this  report 
(see  Figure  6).  In  this  figure,  it  is  important  to  compare  the  results  for  vari¬ 
ous  fuels  at  the  same  temperature.  Thus  one  observes  that  acetylene  is  not  one 
of  the  most  prolific  sooters.  Its  tendency  to  soot  only  appears  great  because 
its  ordinary  flame  temperature  is  so  high. 

IVe  observe  that  the  and  olefins  and  diolefins  have  the  greater  tenden¬ 
cy  to  soot.  The  results  obtained  for  toluene  fall  above  the  olefins  and  indicate 
a  greater  tendency  to  soot.  This  observation  provides  some  insight  into  the  chem¬ 
ical  mechanism  of  soot  formation  and  is  consistent  with  some  of  our  early  propos- 

.  4 

als . 

The  practical  importance  of  these  results  is  that  they  stress  the  great  im¬ 
portance  of  the  stoichiometric  flame  temperature  in  sooting  in  diffusion  flames. 
Recall  there  are  two  temperatures  in  diffusion  flame  combustion  system:  the 
stoichiometric  temperature  of  the  actual  flame  and  the  final  (lower)  temperature 
which  corresponds  to  the  overall  stoichiometry  determined  from  the  fuel  and  air 
flow  rates.  We  shall  address  this  important  point  again  later. 

Numerous  other  experiments  were  performed  to  evaluate  the  effect  of  oxygen 
concentration  on  sooting  diffusion  flames.  The  results  are  shown  in  Figure  7, 
in  which  the  sooting  height  is  plotted  as  a  function  of  the  oxygen  index  (oxygen 
mole  fraction  in  oxidizer  stream) .  The  fact  that  the  tendency  to  soot  undergoes 
a  minimum  (maximum  sooting  height)  is  indicative  of  two  competing  processes  in 
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the  test  procedure.  As  the  oxygen  concentration  rises,  so  does  the  flame  tem¬ 
perature,  pyrolysis  and  the  amount  of  soot  formed.  But  the  rate  of  soot  of 
burn-up  also  increases.  The  two  processes  would  have  two  different  activation 
energies:  high  for  pyrolysis  process  and  low  for  the  burn-up  process.  At  low 
oxygen  concentration  as  the  oxygen  index  is  raised,  the  rate  of  burn-up  dominates 
and  thus  the  flame  height  for  sooting  must  increase.  As  the  concentration  is 
increased  further,  the  flame  temperature  becomes  sufficiently  large  that  the 
high  activation  energy  process  dominates.  At  this  point  the  rate  of  fuel  py¬ 
rolysis,  which  increases  with  temperature,  dominates  and  the  sooting  tendency 
increases  as  the  oxygen  index  increases.  This  model  suggests  that  a  more  ap¬ 
propriate  manner  to  plot  these  data  would  have  been  sooting  height  versus  the 
calculated  adiabatic  flame  temperature  corresponding  to  the  experimental  oxygen 
index.  Different  curves  of  height  versus  index  were  obtained  when  nitrogen  and 
argon  were  used  as  the  inert  in  the  oxidizer  stream.  However  when  the  data  were 
plotted  as  a  function  of  the  temperature  the  nitrogen  and  argon  results  fell  on 

9 

the  same  curve.  The  details  are  given  in  another  recently  published  manuscript. 

These  results  with  respect  to  oxygen  index  may  appear  to  be  in  conflict 
with  other  results  in  the  literature.  Others  found  and  reported  that  the  sooting 
tendency  always  decreased  with  increasing  oxygen  index.  The  anomaly  is  due  to 
the  sooting  criterion  used.  We  used  the  sooting  height  when  the  soot  broke  through 
the  luminous  zone,  whereas  others  used  the  presence  of  soot  at  the  end  of  the 
containing  cylinder  in  which  the  time  for  particle  burn-up  was  extremely  long. 

When  we  used  the  same  criterion,  we  found  the  same  trends  as  others.  However, 
our  results  gave  the  important  practical  result  to  limit  the  amount  of  soot,  one 
must  reduce  the  flame  temperature  and  one  means  for  doing  this  is  to  reduce  the 
oxygen  mole  fraction  by  vitiation.  Again,  the  published  paper  has  more  extensive 
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details . 

During  the  subject  contract  year  we  also  constructed  our  pre-mixed  flame 
experimental  apparatus,  but  sufficient  data  had  not  been  obtained  to  report  any 
results  of  significance. 

IV.  SENSITIVITY  ANALYSIS  STUDIES  OF  COMBUSTION  SYSTEMS* 


Combustion  phenomena  inherently  involve  a  complex  interaction  of  many  sim¬ 
ultaneous  chemical  reactions  and  associated  transport  phenomena.  A  central  ques¬ 
tion  in  this  area  concerns  which  features  of  a  chemical  mechanism  or  transport 
coefficients  are  responsible  for  observable  events  in  the  laboratory  or  on  an 
engineering  scale.  This  general  question  has  been  very  difficult  to  address  due 
to  the  complex  interaction  between  the  system  chemical  and  transport  components. 
In  recognition  of  this  problem,  we  have  been  developing  a  powerful  family  of  sen¬ 
sitivity  analysis  tools  to  quantitatively  address  these  questions.10,11  These 
tools  are  most  useful  in  conjunction  with  on-line  experiments  or  for  the  purpose 
of  designing  such  experiments,  and  the  theory  has  been  tailored  with  this  goal 
in  mind. 

In  its  most  elementary  form  sensitivity  analysis  concerns  an  examination  of 
3c. (x,t) 


the  system  gradients. 


3a. 

J 


where  c^  (x,t)  is  the  i-th  species  concentration 


at  position  x  and  time  t  and  is  the  j -th  system  parameter  (e.g.,  a  rate  con¬ 
stant  or  transport  coefficient).  Other  bulk  thermodynamic  observables  could  also 
be  considered.  An  examination  of  these  gradients  gives  an  immediate  measure  of 
the  relationship  between  the  output  observables  and  the  input  system  parameters. 
More  generally,  an  entire  family  of  gradients  can  be  generated  to  address  a  host 
of  physical  questions.  The  mathematical  structure  of  this  theory  is  very  similar 
to  thermodynamics  and  the  table  below  gives  a  summary  of  the  general  types  of  ques- 
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Contributed  by  Prof.  H.  Rabitz,  Dept,  of  Chemistry,  Princeton  University. 
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tions  that  can  be  addressed. 


4. 


Applications  of  Sensitivity  Analysis 


1 .  Elementary  Sensitivities 


3c^  (x,t) 

9cT 

J 


9a,  (x,t) 

Parameter -Observation  Interdependence  — — - -  .  These  gradients  are 

i 

generally  not  the  inverse  of  those  in  item  1  since  different  dependent  and 
independent  variables  may  be  involved. 


Parameter  Interdependence 


9a- 

i 

9a . 

J 


A  parameter  interdependence  develops 


when  certain  observations  are  assumed  to  be  performed. 

9c.  (x,t ;x; t ' ) 

Interdependence  of  Different  Observations  — 


9Cj 


Functional  Sensitivity  Analysis 


(x,t,x’,t)  6ai  (x,t,x’,t') 


6a 


3 


6a . 

J 


,  etc . 


This  is  the  most  general  case  whereby  the  parameters  themselves  may  be  con¬ 
sidered  functions  of  the  system  coordinates.  All  items  in  the  table  may 
be  extended  to  the  case  of  functional  derivatives. 


32c. 


Higher  Order  Sensitivities 


do. .  da. 

3  k 


,  etc . 


9c.  (x,t) 

Memory  Effects  and  Stability  Analysis  - - — r~rTT  •  This  gradient  is  a 

- = -  oCj  lx  ,t  j 

generalized  memory  function  which  measures  the  stability  of  the  i-th  ob¬ 
servable  at  time  ti  with  respect  to  variations  of  the  j -th  observable  at 
prior  time  t'.  This  function  should  be  distinguished  from  questions  of  ob¬ 
servational  interdependence  addressed  in  ite.n  4  which  assumes  certain  obser¬ 
vations  have  actually  been  made. 

9  JC  ^  t 

Position  and  Time  as  Dependent  Variables  and  i— —  .  Critical  dis- 

3  " "  j 

tances  or  induction  times  can  legitimately  be  considered  functions  of  the 
system  parameters. 

Sensitivity  of  Objective  Functions.  Often  a  desired  physical  observable  is 
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in  fact  a  functional  of  the  output  vector  c_.  Sensitivities  to  the  desired 
objective  function  can  readily  be  generated. 

10.  Sensitivity  to  Missing  Model  Components.  The  sensitivity  of  parameters 
left  out  of  the  model  may  readily  be  calculated  when  evaluating  the  gradi¬ 
ents  at  the  nominal  null  value  of  the  parameters. 

11.  Structural  Sensitivity  Analysis.  Sensitivity  of  global  features  or  struc¬ 
tures  in  the  observables  (e.g.  frequencies  of  oscillating  reactions,  posi¬ 
tions  or  widths  of  chemical  concentration  maxima,  etc.)  can  be  readily  cal¬ 
culated  . 

12.  Mapping  Out  Parameter  Space .  The  sensitivity  gradients  can  be  used  to  gen¬ 
erate  non-linear  interpolation  formulas  for  probing  the  behavior  of  the 
multi -dimensional  parameter  space  governing  a  class  of  chemical/physica 1 
problems . 

1 3 .  Statistical  Error  Analysis. 

The  above  table  represents  the  current  state  of  the  art  in  sensitivity  theory 
and  continuing  research  is  underway  to  both  extend  the  usefulness  of  sensitivity 
theory  as  well  as  carry  out  practical  calculations  for  combustion  oriented  prob¬ 
lems.  In  the  latter  case,  one  particular  area  of  focus  cot  r  r-  a  pro'-x’V  of  the 
interrelationship  between  transport  phenomena  and  chemical  mechanisms  which  is  ad¬ 
dressed  by  item  3  in  the  table.  In  order  to  bring  these  tools  to  a  wider  audience, 

12 

we  are  also  in  the  process  of  developing  a  streamlined  sensitivity  code  which 
would  be  compatible  with  a  variety  of  existing  modelling  packages.  This  work 
should  allow  for  the  addressing  of  a  rich  variety  of  chemical/physical  combustion 
problems  which  were  heretofore  impossible  to  deal  with. 
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